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Abstract 

The Michel parameters of the leptonic r decays are measured using the OPAL detector at LEP. The 
parameters p£, (with £ = e, p) and rj^ are extracted from the energy spectra of the charged 

decay leptons and from their energy-energy correlations. A new method involving a global likelihood 
fit of Monte Carlo generated events with complete detector simulation and background treatment has 
been applied to the data recorded at center-of-mass energies close to y/s = mzo corresponding to an 
integrated luminosity of 155 pb~^. If e-// universality is assumed and inferring the r polarization from 
neutral current data, the measured Michel parameters are: 

p = 0.781 ±0.028 ±0.018, ^ = 0.98 ± 0.22 ± 0.10, 

Tj = 0.027 ±0.055 ±0.005, ^5 = 0.65 ± 0.14 ± 0.07, 

where the value of rj has been constrained using the published OPAL measurements of the leptonic 
branching ratios and the r lifetime. Limits on non-standard coupling constants and on the masses 
of new gauge bosons are obtained. The results are in agreement with the V—A prediction of the 
Standard Model. 
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1 Introduction 


A measurement of the Michel parameters in r decays is presented which involes a novel method to 
fit the energy spectra and energy-energy correlations of the charged decay leptons. The r-pair data 
set used was produced in e'’'e“ collisions at center-of-mass energies close to ^/s = m^o, corresponding 
to an integrated luminosity of 155 pb~^. The parameters are fitted to the lepton spectra of three 
event classes, electron-hadron, muon-hadron and electron-muon, depending on the decay modes of 
the two r decays in the event. Previous measurements of the Michel parameters in r decays exist that 
were performed at LEP and at other e+e“ colliders Unlike previous measurements the 

presented analysis makes use of a binned maximum likelihood ht of fully simulated events to the data. 
It accounts for radiative corrections, detector effects, background processes and selection efficiencies 
in a comparatively original way. 


1.1 Lorentz structure 

In the Standard Model the charged weak interaction is described by the exchange of left-handed W 
bosons, i.e., by a pure vector coupling to only left-handed fermions. Thus, in the low-energy four- 
fermion ansatz, the Lorentz structure of the charged current is predicted to be of the type “V—A ( 8 ) 
V—A”. With this formalism the dominant features of nuclear (5 decay and of ^ decay are correctly 
described |jl^. Deviations from this behavior would indicate new physics and might be caused by 
changes in the W-boson couplings or through interactions mediated by new gauge bosons [^. The 
dominant contribution would then be either the one with the largest coupling or the one being mediated 
by the lightest boson. This means that if there exist contributions to the leptonic decay structure 
other than V—A, for example a right-handed vector coupling arising from a heavy right-handed W 
boson, Wr, or a scalar coupling to a charged Higgs boson, these would emerge first in the decay of the 
massive r lepton. Among all its decay modes, the decays r ^ ei^eW and r —> are the only ones 

in which the electroweak couplings can be probed without disturbance from the strong interaction. 
This makes the purely leptonic r decays an ideal system to study the Lorentz structure of the charged 
weak current. 

The leptonic decay amplitude can be generalized by adding current-current terms for all possible 
bilinear covariants. The most general, derivative-free, four-lepton interaction matrix element for the 
r —> £ decay that is local and Lorentz invariant can be written as (see e.g. 0 ): 

V ^ 7-S,V,T 
e,a;=R,L 

Here 7 denotes the type of the interaction (scalar, vector or tensor) and T^ are 4x4 matrices defined 
in terms of the Dirac matrices: 

rS = l, = = (2) 

The indices u and e denote the chiralities of the r lepton and its charged decay lepton, I", respectively. 
For given (cj, e) the chiralities of the neutrinos are uniquely determined. Tensor interactions exist only 
for opposite chiralities of the charged leptons. This leads to 10 complex coupling constants, for 
which the Standard Model predicts g'LL=I and all others being zero. Choosing the arbitrary phase by 
defining positive leaves 19 real numbers to be determined by experiment. As long 

as one is interested in the relative strengths of the couplings, it is convenient to require the following 
normalization condition: 

^ = i (Ia'llI + I^lrI + \9rl\ + I^rrI ) (3) 

+ (lH^Ll^ + b^Rl^ + b^Ll^ + ll/XRli 

+ 3 (I^lrI + I^rlI ) = 1- 
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This restricts the allowed ranges of the coupling constants to |( 7 ®| < 2, < 1 and The 

overall normalization can be incorporated into Gq which then accounts for deviations from the Fermi 
constant Gp- 


1.2 Michel parameters 

At the Born level, neglecting radiative corrections and terms proportional to only four different 


combinations of these coupling constants, denoted by p, 5 and 77 , determine the shape of the decay 
spectra. In the r rest frame, the leptonic decay width can, for massless neutrinos, be written as: 


d^T 


r— 


dn dx* 


192ir* 

— P-r it COS 9 


{l-x*) + 5t(-x* -2 


mt (1 — X*) 


(4) 


7T * 771^ 

Here x* = ^max is the scaled energy of the charged decay lepton in the r rest frame with ^ 

being its maximal energy and cos 0 * is the angle between the r-spin direction and the momentum of 
the decay lepton. is the average r polarization. After integration over cos 6* and boosting into the 
rest frame, the spectrum has the form; 


H{x) = f{x) + Pr g{x), with 
/(x) = a(x) + pb{x) +r]e{x), g{x) = i c(x) +i6d(x), 


(5) 

( 6 ) 


where /(x) and g(x) describe the isotropic and the r-spin-dependent part, respectively, and a(x)... e(x) 
are known third-order polynomials in the scaled energy x = ^. Hence, allowing the most general 
couplings, the shape of the spectra can be described by the four Michel parameters [IT, ]T| for which 
the Standard Model predicts the values /9 = |,| = 1,(5 = | and g = 0 according to a V—A structure 
of the charged weak current. Their definitions read in detail^ 

P = ^I^llI +^lfl'RRl +y^|5'llI +y^|5'lrI +y^|5'rlI +Y^lfl'RRl 

+ ^|5lrI + ^IffRul “ ^^'^(SlrSlr) “ ^I4e(fl'RLfl'RL) ) 

i = Is'llI +3|5'lrI “3|fi(^L| - Is'rrI +5|5'lrI “ 


+ ^|5llI “ ^Is'lrI + ^I^rlI “ ^Is'rrI + 4Re(gLR9LR) “ 4Re(5'pL5'RL)) 

“^I^lrI +^|9rlI + ;^I^r(9lr9lr) - ^Re(9|L9RD. 


= -Re 


^LL^RR + S'RRS'LL + 5'Rl(5'LR + S^Lr) + fi'LR(5RL + ^^Rl) 




..V / S* 


,T* ' 








(7) 


One of the parameters, g, can also lead to a change in the leptonic decay width of the r lepton. 

If contributions from other couplings exist, they are not necessarily the same for r —> 
and r —> Therefore, the Michel parameters for decays into electron and muon are measured 

independently. Because the ry^-term is suppressed by a factor there is almost no sensitivity to 

^ Since by definition 5 contains a factor i it is convenient to use (^5) instead of 5. 
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r]e from the r ^ decay spectrum. Thus in the following rye is set to zero. This leaves the 7 

parameters pe, Ce, and to be determined. 

To test the Standard Model prediction the parameters will also be fitted under the assumption of 
e-p universality. It is possible to test explicit extensions, for example by focusing on vector couplings, 
or by allowing only one additional scalar contribution. In this way, mass limits for a right-handed W 
boson, Wr, as well as for a charged Higgs boson are determined. 


1.3 Lepton-lepton correlations 

At LEP, r pairs are produced with almost perfect spin correlation.^ This allows one to measure the 
spin-dependent part of the decay spectra with high sensitivity by employing the correlations between 
both r decays in the event 0,0- For parallel r spins the correlation function can be written as: 

I(xi,X 2 ) = f{xi)f{x 2 ) +g{xi)g{x 2 ) - Pr [f{xi)g{x 2 ) + f{x 2 )g{xi)] , (8) 

where f{x) and g{x) are the above third-order polynomials and Pr is the average polarization of 
the T~ lepton. In contrast to uncorrelated single decay spectra (eq. P) here the product of the two 
spin-dependent parts, g{xi)g{x 2 ), without the suppression by the r-polarization. Thus, in case of two 
leptonic decays in the event, the correlation function provides high sensitivity to the parameters ^ and 
^6 whereas the more frequent events with single leptonic decays contribute with high statistics to the 
measurement of p and rj. 


1.4 Tau polarization 


In addition to the charged current couplings which determine the decay of the r lepton, the couplings 
to the neutral current which are responsible for r production influence the shape of the spectra. The 
difference in the left- and right-handed couplings causes the r leptons to be produced polarized 
affecting the spin-dependent part of the decay spectra. In principle, the average r-polarization, Pr, 
could be measured along with the Michel parameters. However, this introduces additional correlations 
between the fit parameters and limits the accuracy of the Michel parameter measurements while it 
does not reveal anything new about the decay structure. Since extensions to the charged sector of the 
weak interaction do not a priori change the neutral current, a different approach is pursued in this 
analysis. Instead of testing both sectors simultaneously at the expense of less sensitivity, the charged 
current is investigated in the most general way, while the neutral current (i.e., the r-polarization) is 
assumed to be described by the Standard Model couplings with adequate accuracy. 

Since all direct measurements of Pr from e"’“e“ ^ data at ^/s = m^o implicitly have assumed 

a V—A coupling in the r decay, the use of these results as input to this analysis would introduce a bias. 
Therefore, measurements of the neutral current which are independent of the charged sector have been 


used to calculate Pr using the zfitter [19| package. The uncertainty arising from this procedure has 
been studied with the systematic errors. Its impact on the Michel parameter measurement is smaller 
than the one that would have been introduced through correlations between the parameters and the 
polarization if Pr were fitted. 

As input to the Pr calculation the following values were used the preliminary LEP measure¬ 
ment of the mass m^o = (91.1884 ± 0.0022) GeV, the CDF/D0 combined value for the top mass 
rrit = (175 ± 6) GeV, mjjo = 300GeV with 60 GeV < mjjo < ITeV, as{mio) = 0.118 ± 0.003 and 
a{mio)~^ = 128.90 ±0.09. These values yield the predicted value of the r-polarization at the peak 
as Pr = —0.1391 , where the dominant uncertainty is due to the unknown Higgs boson mass. In 

this analysis, data taken at ^/s = m^o are used together with data below and above m^o. Averaging 
over the energy dependence of Pr according to the data set used yields the same quoted value. Over 
the full range this leads to an uncertainty in Pr of approximately 5 % which is still smaller than the 
error of the current measurements (see e.g. pl| ). 

^For V and A type couplings in the production, the and r~ have opposite chiralities. 
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2 Event selection 


2.1 OPAL detector 


The OPAL detector is described in detail elsewhere |^, ^|. Here, only a brief summary of the main 
components shall be given. The innermost subdetector is a micro-vertex detector with two layers of 
double-sided silicon-strips. It is enclosed by a system of three different drift chambers, a precision 
vertex chamber with axial and stereo wire readout, a large cylindrical drift volume (jet chamber) 
with 24 azimuthal sectors of 159 signal wires each, and a surrounding set of z-chambers with wires 
perpendicular to the beam direction. The central tracking system is contained inside a magnetic 
coil which provides a solenoidal field of 0.435 Tesla. This leads to a resolution of the transverse 
momentum of (J-ptlVt ~ 1-5 • 10“^ • pt (GeV). In addition, the specific energy loss, dE/dx, of charged 
particles is measured in the jet chamber. The electromagnetic calorimeter (ECAL) is built of 11,704 
lead glass blocks of approximately 25 Xq radiation lengths, providing an energy resolution of typically 
oeIE k, 12 %/-\/E (GeV). In front of the EGAL, a thin gas detector (the presampler) measures 
electromagnetic showers beginning in the material of the inner detector components. The ECAL is 
surrounded by the iron return yoke of the magnet which is instrumented with limited streamer tubes to 
serve as a hadron calorimeter (HCAL). The whole detector is enclosed by 4 layers of muon chambers, 
giving the position and the direction of the penetrating particles. 

As observables for the fit, different variables in the two leptonic decay channels are chosen. In case 
of the r —> decays, the track momentum of the muon measured in the jet chamber, ptrack; is 

taken, whereas for the r —> e decays the energy deposited by the electron in the electromagnetic 
calorimeter, Eduster, is used. The cluster energy includes a large fraction of final state photon radiation 
from the electron and is thus less dependent on the modeling of the photon radiation in the generator 
and detector simulation. Both variables are scaled to the beam energy which is used as the estimator 
for the maximal energy of the decay lepton. 


2.2 Tau pair selection 


Erom the data collected with the OPAL detector during the years 1990-1995, r-pair events are selected 


in several steps. The selection follows the strategy described in earlier OPAL publications |^, ^| and 
given in detail in |^]. Eirst, lepton pairs are preselected by requiring exactly two charged jet-cones 
of 35° half-opening angle with low track and cluster multiplicity. Two-photon events are rejected 
by requiring either a large visible energy or an unbalanced transverse momentum sum, and a small 
acollinearity. 

Erom this sample Bhabha events are removed based on a large sum of cluster energies or a large sum 
of track momenta in conjunction with large cluster energies. Afterwards, /r-pair events are eliminated 
if consistent with high track momenta, small energy deposit in the calorimeter or signals in the muon 
chambers. The remaining 147 042 events are almost entirely r pairs (see section |2.5| for the remaining 
background in the used event classes). The geometrical acceptance of this selection covers the region 
|cos(9| < 0.95. 


2.3 Tau decay mode identification 

A likelihood selection is used to identify the r-decay modes in the two cones. It distinguishes between 
the 1-prong decays r —> eUeU-r, r —> and r —> hir^- where h is either vr/K, /) or ai —> 7r27r^. It 

makes use of a set of variables which allow the discrimination of different channels These variables 
include; 

Eciuster/ptrack “ the ratio of the energy deposited in the electromagnetic calorimeter (EGAL) with 
respect to the track momentum measured in the central detector, 

^Kaons are not separated from pions and are counted among the corresponding 7r-channels in the following. 
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dE/dx - the specific energy loss in the jet chamber, 

-^biolL ” number of ECAL blocks that contain 90 % of the measured energy in the cluster 
associated to the track, 


-^neutral / Eq, 


the fraction of ECAL energy in the cone that is not associated to the track, 


(A(^)max 


the maximum angle between the track and a presampler cluster assigned to the 
cone. 


khpres “ the width of the largest presampler cluster. 

Ambits/-^layers “ the average number of hits per active layer in the hadronic calorimeter (HCAL), 

A^ 7 hits “ the number of hits in the last 3 HCAL layers and in the 4 layers of the muon 
chambers. 


khmuon “ the matching probability between the extrapolated track and a muon chamber track 
segment. 


The measured variables are then compared to a set of reference distributions that have been produced 
by a Monte Carlo simulation of the considered decay modes (see section . 

Based on the variable x*, the expected fraction of decay modes of the type j is given as: 






(9) 


where // are the normalized probability densities taken from the reference distribution for the respec¬ 
tive variable i and Amodes = 5 is the number of considered decay modes. The information from all the 
variables is then combined into the product likelihood for the hypothesis j: 


A^variables 

i=l 


( 10 ) 


Normalized to all considered alternatives the expression yields the relative likelihood of the decay 
being of type j: 


pj = 


a 


^^^modes pj 


( 11 ) 


which in the case of uncorrelated variables can be interpreted as a probability. To obtain high purity 
samples, the cones involving the lepton are required to be identified with a relative likelihood of 
> 90%. The simulated reference distributions are checked against the data using a tagging 
technique. Eor that purpose, the likelihood variables of each detector component are compared between 
data and Monte Carlo based on a clear decay mode identification from the other components. A 
comparison of some of the variables can be found in [^|. 


2.4 Selection efficiency 

Throughout the analysis, the simulated events and the data are treated in an identical manner. Thus, 
the efficiencies for the event selection and the decay mode identification are accounted for in the fit as 
modeled in the Monte Carlo. The validity of the simulation has been tested by comparing Monte Carlo 
efficiencies as a function of energy / momentum with data control samples (see section P). Table I 
lists the efficiencies for the decay mode identification within the fiducial volume of | cos 6\ < 0.95. The 
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Generated mode 

T ^ e 

Identified mode 

T ^ T ^ TTIffr T 

^ pVr 

T — > ai z/,- 

r ^ e UeUr 

96.5 

- 

1.1 

1.7 

0.7 

T p.VpVr 

0.3 

92.2 

7.0 

0.5 

- 

T ^ TTUt- 

1.4 

3.3 

82.1 

9.5 

3.6 

T ^ pVr 

0.3 

0.2 

10.3 

67.8 

21.3 

r — > ai 

0.1 

- 

1.3 

25.8 

72.7 


Table 1: Efficiencies in % for the likelihood identification of 1-prong r decay modes within the fiducial 
region. The diagonal values are the efficiencies for the true identified modes whereas the off-diagonal 
values represent the misidentification probabilities. 


inefficiency of the muon channel originates mainly from the low x region where the separation from 
pions becomes poorer. 

Based on the decays of both r leptons, the events are divided into mutually exclusive samples: 
lepton-lepton correlations (e-/r) and single-lepton decays (e-hadrons, /r-hadrons). As mentioned in 
section the correlation spectra have a high sensitivity to the spin-dependent parameters (^, 
while the single decays provide high statistics for the isotropic parameters (p, rf). 

Decay correlations in e-e and p-p events are subject to large background contamination from 
e+e“ ^ e+e” and e+e" —> pffpi~ processes as well as from two-photon events 77 ^ e+e" and 
77 —> These background sources distort the correlation spectra in the two most sensitive 

regions where both leptons have high energies or both have low energies, respectively. Due to the 
uncertainties in the estimation of these backgrounds, the e-e and p-p event classes would contribute 
with a dominant systematic error to the fit result. For this reason they are not included in this 
analysis. 

In table |2|the overall efficiencies for the three event classes are given. The additional inefficiencies 
with respect to the decay mode identification arise from geometrical cuts that are performed on the 
individual event classes to exclude insensitive or inadequately simulated regions of the detector. These 
regions include the proximity of the anode wire planes of the jet chamber, (/>sector < 0.4°, the extreme 
forward region of the electromagnetic calorimeter, | cos0| > 0.9, and a few small regions not covered 
by the muon chambers. 




Event Class 



Tl ^ eVeVr 

Tl ^ P UpUr 

n ^ eUeUr 


T2 —> h Ur 

T2 Ur 

T2^ P 

Sample size 

19369 

21190 

5834 

Efficiency 

83.1% 

88.6% 

85.8% 


Table 2 : Selection efficiencies for the three considered event classes. 

Figures |l|(a) and ||(b) show the selection efficiencies for the single leptonic decay samples e-h and 
/i-h as a function of the scaled energy xe and scaled momentum Xp, respectively. 

2.5 Background 

Various sources of background in the r-pair sample and in the r —> and r —> samples 

have been investigated. The main contribution comes from misidentification of the hadronic r decays. 
Due to the distinctive signatures of electrons and muons in the detector, cross-talk between the two 
channels is negligible. The main backgrounds for the r —> ei^e^^r mode originate from r —> ttz^t and 
T —> decays. For the r —> mode the dominant background source are r —> ttz/t- decays. 






Table ^ lists the different backgrounds in the three event samples. For simplicity, the three identified 
hadron channels have been summed up. For the single decay spectra, double lepton events in which 
the recoil lepton has been misidentified as a hadron are also quoted as background. 


Background source 

Event class 
e-h //-h e-fj. 

leptonic r e-e 

e-pL 

/i-e 

n-n 

0.73% - 0.00% 

1.70% 

0.99% 

1.64% 0.12% 

hadronic r h-h 

e-h 

h-fi 

0.72% 1.55% 

1.73% 

0.78% 

non-r e e —> e e 

ee —> /r/i 

77 ^ ee 

77 ^ ///i 

0.06% - 0.00% 

0.56% 0.06% 

0.03% - 0.00% 

0.21% 0.04% 

Total 

3.3% 5.0% 2.8% 


Table 3: Identified events and background from r decays, Bhabha, ^-pair and two-photon events. For 
simplicity the three hadronic channels have been summed up. 

Furthermore, Bhabha and /r-pair events may pass the selection and contaminate the samples. This 
can occur, for example, for an e^e“ —> e'’'e“ event where one electron is misidentified as r —> 
decay, or to an e'’'e“ ^ event where one muon fakes a r —> vr z/,- decay. 

As mentioned earlier, the presented method allows to account for all backgrounds that are simu¬ 
lated in the Monte Carlo. Since either the background shapes do not depend on the Michel parameters 
(hadronic and non-r background) or their dependence has, due to the small fraction, no significant 
impact on the shape of the signal distribution, the residual events can be independently added to the 
fit spectrum. 

3 Fitting method 

To extract the Michel parameters from the observed spectra, a binned maximum likelihood fit of a 
set of Monte Carlo spectra has been applied to the data. Compared to fits that involve analytical 
functions to describe the data distributions, this method has several advantages. It includes radiative 
corrections at the generator level where their description is more easily accessible than it is through 
the convolution of photon radiation probability functions. It provides a full simulation of the detector 
response and thus accounts for the energy and momentum resolution as well as for the selection 
efficiency in an elegant way. It also accounts for contamination from misidentified r decays and from 
other background sources, like Bhabha events, /i-pair events and two-photon events. This makes it 
unnecessary to unfold resolution and efficiency effects from the data. A particular benefit is that 
selection criteria which restrict the phase space do not present a problem for the fit because all 
requirements can be placed identically on data and Monte Carlo. In addition, observables for each 
decay mode can be chosen to accommodate the particular capabilities of the detector rather than to 
facilitate the theoretical description. 
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3.1 Linear combination 

Since the single decay spectra depend linearly on the Michel parameters, it is possible to decompose 
any observed spectrum into different basis spectra with each representing a specific set of parameters. 
To describe the whole parameter space including the constant term, five spectra have to be mixed with 
coefficients that add up to unity. For a particular choice of the basis parameter sets, the corresponding 
coefficients can be calculated by solving the respective equation system. This allows one to determine 
the values of the 4 Michel parameters by fitting the coefficients with respect to a given basis (i.e. the 
relative contributions of the basis spectra) to the data. 

In the case of the energy-energy correlation spectra between two leptonic decays in each event 
(double decay spectra) this method is still applicable. These spectra can be represented by a compo¬ 
sition that is bilinear in the 2x4 Michel parameters of the r decays into electron and muon (or of 
second order in the 4 parameters if the leptons in both hemispheres are identical). The correlation 
basis is then the tensor product of two single decay bases. Now, 25 coefficients appear in the decom¬ 
position, but are not all independent. From these, the 2x4 Michel parameters for e and /r decays, 
respectively, can again be calculated by matrix inversion (see section |3.2| ). In this formalism, the e-e 
and /i-// correlations with only 4 free Michel parameters are just a special case of the e-// correlations .0 
They will not be discussed explicitly since they are not used for the analysis. 

With this method, which is based on Monte Carlo event generation, it is possible to describe any 
value of the Michel parameters with a finite sample of events by varying the appropriate contribution 
to the spectrum. 


3.2 Parameter basis 


To describe the entire Michel parameter space, a basis has to be chosen which accounts for any possible 
combination. It is obvious that the parameter sets corresponding to pure S,P,V,A,T couplings do not 
fulfil this requirement because they do not involve any interference terms between different couplings. 
Also, the canonical basis with (p, ^6, rj) = ( 1 , 0,0, 0), (0,1,0, 0),... etc. is not physically meaningful, 

because it is not possible to have ^5=1 and all other Michel parameters zero. However, it is desirable to 
have a “nearly orthogonal” set in order to avoid large correlations between the coefficients. Examples 


of different couplings and the corresponding Michel parameters can be found in table 1C. 
Taking this into consideration, the following basis has been chosen: 


/ P \ 
\ V J 


= Cl 


f 3/4 \ 

1 

3/4 

V 0 y 


+ C 2 


/ 1 \ 
0 
0 

voy 


+ C3 


/ 0 \ 
1 
0 

voy 


+C 4 


/ 3/4 \ 

0 

0 

Vi/2y 


+C 5 


/ 0 \ 
0 
0 

voy 


( 12 ) 


with = 1- Here, each vector on the right-hand side stands for a spectrum generated with the 

quoted set of Michel parameters. 

In a more general notation^, one can write for the case of the single lepton spectra: 


cie = Me- C£, 


(13) 


where is the vector (1, pe, ^e, 6e, rje), ce is the vector of coefficients (ci, C 2 , C 3 , C 4 , C 5 ), and Me is a 5 x 5 
matrix given by 


Me = 


/ 

1 

1 

1 

1 

1 

\ 


Pll 

Pe2 

Pes 

Pu 

Pn 





in 

iu 

in 



5 11 

<5f2 

^n 


^n 


V 

Vei 

Pi2 

Pn 

PU 

Pn 

y 


( 14 ) 


“^They are formally analog to e-ii-correlations under the assumption of e-^i universality. 
®To simplify the notation, 5 is written instead of ^5 in this section. 
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Here pii denotes the value of in the first basis spectrum, p £2 its value in the second spectrum and 
so forth. Hence, for a given set of Michel parameters the coefficients with respect to the basis can 
be calculated from 

C£ = (15) 

Analogously, for the case of the double e-p spectra one can write: 

• Ce^, ( 16 ) 

where the vector qe^ is the outer product of the vectors qe and q^ ( gef^j = Qei ' g^lj )> is the 
corresponding outer product of the matrices Mg and M^, and Cg^ = (ci,...,C 25 ) is the coefficient 
vector. Again, for a given set of 8 Michel parameters the corresponding 25 coefficients in the actual 
basis are obtained by multiplying the vector qg^ = (1, pg, ^g, 5g, pg, p^, PeP^, ^gP^,..., Pgp^) with the 
inverse of the 25 x 25 matrix Afg^. 


3.3 Monte Carlo simulated samples 


The r-pair Monte Carlo sample was generated using the koralz-3.8 |27] generator and a modified 
version of the tauola-1.5 |^| decay library which was extended to include the full generalized matrix 
element This allows the variation of the Michel parameters over their whole range. 

The tauola-1.5 version was preferred over the newer version tauola-2.4 |^,M| since it has 
a more general approach to real photonic corrections. The tauola-1.5 library uses a factorization 
ansatz in order to produce real photons, where each charged particle radiates independently of its 
decay matrix element by internally applying the photos package. On the other hand, the newer 
version TAUOLA-2.4 contains full 0{a) corrections to take the interference between photons radiated 
from the r and photons radiated from the decay lepton into account. Although this might be more 
precise in the low energy regime, it assumes a V—A type of interaction and is thus model dependent. 
The former version, TAUOLA-1.5, is capable of producing photons for all possible decay structures at 
a reasonable level of accuracy. It has therefore been chosen for this analysis. 

Since finite Monte Carlo samples are used to describe possibly small variations in the shape of the 
spectra, it is vital to keep statistical fluctuations under control. If the distinct spectra were generated 
independently their differences would be smeared by Gaussian errors. In particular, it could happen 
that in a certain bin, the theoretical prediction for one spectrum is higher than for another but the 
generated number of events is lower. To avoid such fluctuations, each event is used for as many spectra 
as possible. This means that most of the Monte Carlo spectra share a large fraction of common events. 
An acceptance/rejection method is used where an event is flagged as accepted for each spectrum for 
which the generated random weight is below the prediction (and not just the standard V—A) and it 
is rejected only if it belongs to none of the considered spectra. This procedure guarantees that the 
difference of any two spectra has the right sign in all bins. As a side effect, it makes the generation of 
the Monte Carlo samples much more efficient. 

Non-r background sources such as Bhabha, /r-pair and two-photon events have been generated 
using the babamc 


generators, respectively. The re- 


KORALZ Q and VERMASEREN 


sponse of the detector to generated particles is modeled using a simulation program |]^ based on the 
GEANT package. The simulated detector response has been checked with various control samples 
(see section El- 


3.4 Fitting procedure 

The appropriate Monte Carlo distributions, each representing a basis parameter set, are prepared 
with the identical binning as used for the observed lepton spectra. The coefficients of the spectra, 
C£, are determined using a binned maximum likelihood fit to the data. To avoid any dependence on 
the description of the overall efficiency, no constraint is made on the overall normalization. Since 
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the coefficients are not independent, they are not varied themselves, but are calculated from the 
corresponding Michel parameters, as has been explained above. Thus, the Michel parameters are 
varied and fitted directly. The actual minimization and the determination of the covariance matrix is 
performed using the minuit package. In doing this, a likelihood is computed for every mixture, 
assuming Poisson errors in each bin. Although the generated Monte Carlo sample is roughly four 
times larger than the data sample, some spectra still have bins with only few entries. This occurs 
in particular in the correlation distributions in regions of high {xi,X 2 ) as the various couplings may 
differ drastically for extreme momentum configurations. It is known that, in the case of small bin 
entries, ignoring the Monte Carlo errors biases the mean value of the fit and underestimates its spread. 
Therefore, fluctuations of both data and Monte Carlo have been taken into account. To accomplish 
this, an adjusted likelihood is calculated |M| by finding in each bin i the most probable expectation 
with which data and Monte Carlo are consistent: 


.^bins .^spectra .^bins 

ln£ = ^ ^ {aji \nAji - Aji) + ^ {di In fi - fi). (17) 

i=l j=l 2=1 

Here Agpectra is the number of Monte Carlo basis spectra (5 or 25), di is the observed number of data 
events, aji is the generated number of Monte Carlo events in spectrum j and Aji is the best estimator 
for the Monte Carlo in the light of the data. The Monte Carlo expectation fi = Y^jPjAji is the 
composition of the best estimators using the mixing coefficients pj. The first term of the adjusted 
likelihood (^) accounts for the agreement between the actual Monte Carlo distribution {aji) and the 
ideal distribution {Aji) for all spectra j. The second term accounts for the agreement between the ideal 
composition (/j) and the data {di). The likelihood is maximized with respect to both the coefficients 
Pj and the estimators Aji.^ In the case of the double lepton spectrum, the index i is replaced by 
two indices. For clarity, all the following expressions correspond to the single lepton spectra. The 
generalization to the double lepton case is straightforward. 

The above equation represents the correct treatment of the problem, provided that the generated 
numbers, aji, are statistically independent. As described before, this is not the case for the prepared 
Monte Carlo spectra because they have (most) events in common. It is, however, possible to rewrite 
the expectation in each individual bin by means of independent numbers. To this end, the spectra are 
ordered by increasing numbers of events in the considered bin, and the coefficients are recalculated in 
terms of the differences between the bins as follows. If one writes for a particular composition (setting 
boi ■= 0 ): 

-^spectra -^spectra 

fi = ^ Cj bji = ^ Pj {bji — bj-ij), (18) 

i=i i=i 

where bji is the generated number of events in bin i of basis spectrum j, then the differences aji = 
bji — bj-ij are all independent. It follows that: 


-^spectra 

Pj = ^ Ck, 

k=j 


(19) 


where the coefficients are expressed in terms of the Michel parameters as described above (sec¬ 
tion 3.2): 

Ck = {M~^)ki ■ Qi- ( 20 ) 

After this transformation the above expression for the likelihood is applied. 

®For each set of coefhcients Pj the numbers Aji can be calculated by solving an equation system. 
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3.5 Constraint on the eta parameter 

It has been mentioned in section |1.2| that the Michel parameter r] corresponds to a change in the 
partial decay width when it deviates from zero. However, to eliminate any dependence on the overall 
efficiency, Monte Carlo predictions are always normalized to the data. Due to the large correlation of 
rj with the p parameter, configurations are possible where the spectral shape is changed only slightly 
while the value of p is inconsistent with the observed branching ratio.[^ 

Such a situation is avoided by constraining the branching ratios to their measured values through¬ 
out the fit. In the general ansatz for the Lorentz structure, the leptonic width is in lowest order 
changed to: 

T,{p,) = rf (^1 + , (21) 

where is the Standard Model width (see e.g. |^). From the measurement of the r lifetime, r^, 

the expected branching ratio, = B{t ii/iUr), depending on the value of pi can be calculated as: 

-B£(%) = r^(%)r^. (22) 


This relation can be used to calculate the most probable value of p, 

^ 1 rrir I Bi 


m = 


4 mi 




-1 


Using the published OPAL results for the r —> branching ratio, B^ = 0.1736 ± 0.0027 


the T lifetime, Tr = (289.2 ± 2.1) fs |^], one determines: 

p-^ = 0.032 ± 0.073, 


(23) 
and 

(24) 


which is consistent with zero. The constraint to p is applied by adding the following term to the log 
likelihood: 


1 rconstraint _ B() 

' “ 2 (r,(r,,)Ar.)2 + (AH,)2- 


(25) 


The use of this constraint makes the fit result for dominated by the branching ratio and lifetime 
measurements while the other parameters are still sensitive to the allowed variation in the ? 7 ^-dependent 
part of the shape. 


3.6 Checks of the fitting method 

To verify the reliability of the fit method, various checks have been performed using the Monte Carlo 
event samples. First, it has been tested that the fit of the linear composition can reproduce the genuine 
parameters of a specifically generated sample. It has also been proven that any genuine spectrum and 
its corresponding mixture of basis spectra are consistent within the statistical errors. 

Second, the statistical errors of the fit parameters have been checked. For that purpose, the 
Monte Carlo sample has been divided into several pairs of subsamples with each pair representing a fit 
sample and a fake data sample. Then the fit has been performed for each subset separately, and the 
distribution of the fit result has been checked. It has been verified that, when applying the adjusted 
likelihood, the individual errors are consistent with the spread of the mean values. 

Third, it has been checked whether the particular choice of the parameter basis affects the fit 
results. Other sets of basis spectra than the quoted one have been used to describe a given Monte 
Carlo distribution. While some bases turned out to be less sensitive to one or more of the parameters, 
all alternatives have been found to be consistent within their fit errors. 

^ It has been shown that the leptonic branching ratios provide a sensitive observable on the Michel parameter p . 
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4 Results 


The result of the most general global fit to single-e, single-// and e-/i-correlations with seven free 
parameters is shown in table Table gives the individual correlation coefficients. The first 



r —> e Ugiyr 


T ^ p 

Pe 

0.779 ±0.047 ±0.029 

Pu 

0.777 ± 0.044 ± 0.016 


I.I3± 0.39 ±0.14 


0.79 ±0.41 ±0.05 

(e5)e 

0.72 ±0.31 ±0.14 

(e<5)p 

0.63 ±0.23 ±0.09 

Ve 

0 (fixed) 


0.010 ±0.065 ±0.001 


Table 4: Results of the global fit to the e-h, /i-h and e-/x energy spectra. The first error reflects both 
data and Monte Carlo statistics, the second error is systematic (see section fA). The error on 
includes the errors on the muonic branching ratio and the r lifetime. 




me 

Pu 



Vu 

Pe 

-0.716 

-0.700 

-0.249 

0.229 

0.344 

-0.035 



0.592 

0.251 

-0.337 

-0.307 

-0.006 

me 



0.321 

-0.289 

-0.442 

0.037 

Pu 




-0.427 

-0.590 

0.461 

Cu 





0.271 

0.132 

mu 






-0.116 


Table 5: Correlation coefficients between the parameters of the global fit. 

errors are due to data and Monte Carlo statistics as obtained from the fit, the second ones are due to 
systematic uncertainties (see section ^T|) . Figures and |3| show the single decay spectra for r —> e 
and T ^ fi decays, respectively. In both figures, the light shaded histogram is the adjusted Monte 
Carlo prediction using the Michel parameters from the global fit (table 

Note that the fit range for the r —> e spectrum extends from x = 0 to 0.9 while for r —> /t 
as well as for the e-//-correlation spectra the full range (x = 0 to 1.0) has been used. This is, because 
the background from Bhabha events causes a relatively large systematic uncertainty at high x to single 
r —> ei^ei^T decays which is no longer present in the correlation spectra where the recoil is required to 
be a r —> // decay. For the r —> // single decay (fig. the e+e" —> pfi contamination at 

high X is well simulated and no restriction of the fit range is necessary. 

In the plots below the single lepton spectra the differences between data and Monte Carlo are 
shown. Since only the shape of the likelihood function is significant and not its scale, there is no 
absolute criterion for the quality of the fit from the maximum likelihood method itself. As an estimator 
for the confidence level, the probability^ for the plotted bins is appended to figures ^ and 

Figures ^ and |5| are two different representations of the e-//-correlation spectrum, displaying slices 
of Xe for fixed x^ (fig. ^ and vice versa (fig. ^). In all figures the dark shaded histogram represents 
the total background from misidentified r decays and the black one the corresponding fraction of 
e+e“ —> e'''e“ and e^e“ —> p~ events. The background decomposition is listed in table ^ The 

® It has to be noted that this yfi probability is biased towards higher values. Since the number of fit parameters that 
constrain a specific spectrum cannot be assigned to the different event classes in a unique way the (number of bins - 1) 
has been taken as the number of degrees of freedom. 
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fitted Michel parameters of table 0 are consistent with the V—A prediction of the Standard Model 
and describe the spectra in figures |2[j^ well. 

4.1 Systematic errors 

Various systematic uncertainties arising from different sources, both detector specific and inherent to 
the method, have been investigated. Besides these systematics, uncertainties concerning the Monte 
Carlo simulation at the generator level, i.e., radiation effects, branching fractions, process kinematics 
etc., have been found to be negligible. The Monte Carlo statistical error is already reflected by the fit 
errors as explained above. 

For the study of systematic uncertainties imposed by the Monte Carlo simulation various data 
control samples have been prepared. Two samples, one of Bhabha and one of /r-pair events, have 
been selected using the criteria mentioned in section |2.2| . The energy and momentum distributions 
of these lepton-pair samples show narrow peaks at the beam energy. Two other samples with single 
lepton cones have been prepared using a tagging technique. This has been done with preselected two- 
photon events, 77 ^ e+e" and 77 —> where one lepton is tagged based on a clear decay mode 

identification so that the other cone can be investigated. The same has been done with e'’'e“ —> e'’'e“ 
and e+e" —> events. 

The following systematic effects, summarized in table P, remain: 

- The absolute scale of the energy and momentum measurement of electrons and muons, respec¬ 

tively, has been varied in the Monte Carlo. The scale uncertainty has been determined by 
comparing e'*“e“ —> e'*“e“ and e"*“e“ —> events from the data control samples with the 

corresponding Monte Carlo events. Rescaling the energies and momenta of the Monte Carlo 
events used in the fit within the observed uncertainty of 0.3% and 0.2%, respectively, results in 
changes of the Michel parameters as given in table 

- The uncertaintiy in the energy and momentum resolution has been determined from the same 
control samples which have been used for the previous study. Variation of the resolution leads 
to the quoted changes. 

- The efficiency as determined from the Monte Carlo has been compared to the tagged control 
samples that include Bhabha and //-pair events covering the high x region as well as two-photon 
events at low x. The events have been used by tagging one lepton while examining the efficiency 
for identifying the other lepton. The ratio of the efficiencies in data and Monte Carlo has been 
found to be consistent with unity. The uncertainty on this efficiency ratio has been estimated 
by fitting the low and the high x values with a straight line and determining the la bounds of 
a possible slope. Since the overall efficiency is not used in the fit only an energy/momentum- 
dependent discrepancy could possibly affect the result. This has been accounted for by weighting 
the Monte Carlo events within the uncertainty of the possible slopes. 

- The background contributions from r and non-r sources which are estimated from the Monte 
Carlo simulation have been varied in the fit. Previous studies of the decay mode identification 
mentioned in section |2.3| showed that the reference distributions of the likelihood variables agree 
with those of tagged data samples fairly well. Reweighting the Monte Carlo such that it per¬ 
fectly matches the data causes only slight changes in the resulting purities. The uncertainty of 
the background is savely estimated to be at the 10 % level, with the exception of the Bhabha 
background which has been varied by a factor of 2. The reason for this is a relatively poor 
modeling of the Bhabha background especially in the forward region. Due to the small fraction 
of the background, the variation had no significant impact on the fit result. 

- The dependence of the fit result on the energy and momentum range has been investigated 
by omitting the outer bins at low and high x values. For the single // and the e-fi spectra the 
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variation of the fit with and without the high x bins also reflects the sensitivity to the description 
of the e"’“e“ —> /i"*"//” background. No statistically significant change has been observed. 

- Different sets of basis spectra have been used to study the influence of the particular choice of 
the Michel parameter basis on the description of the data. The same behavior has been observed 
as for the Monte Carlo study mentioned above with the cross checks of the fitting method. Due 
to the finite Monte Carlo statistics some bases were less sensitive to the Michel parameters than 
the preferred one, however, all alternative fits were consistent within their respective fit errors. 
This means that the basis choice does not contribute to the systematic error. 



Ape 

<1 

A(e5)e 


<1 

A(e5)^ 

Ar?^ 

Energy scale 

0.017 

0.11 

0.11 

- 

- 

- 

- 

Momentum scale 

- 

- 

- 

0.013 

0.04 

0.04 

0.001 

Energy resolution 

0.003 

0.03 

0.04 

- 

- 

- 

- 

Momentum resolution 

— 

— 

— 

0.004 

0.05 

0.02 

0.001 

Energy-dependent efficiency 

0.023 

0.08 

0.08 

- 

- 

- 

- 

Momentum-dependent efficiency 

- 

- 

- 

0.009 

0.07 

0.01 

- 

Total 

0.029 

0.14 

0.14 

0.016 

0.09 

0.05 

0.001 


Table 6 : Contributions of systematic uncertainties to the error of each Michel parameter for the global 
fit. A hyphen indicates that the listed effect contributes less than 0.001 (p,ri) or 0.01 to the 

Michel parameter error. Note that the value of is dominated by the branching ratio constraint. 


4.2 Lepton universality 

With the assumption of universality between electron and muon, i.e., with all couplings gf^ being 
the same for r —> et'eW and r —> one set of Michel parameters can be used to describe both 

leptonic decays. The fit then yields the results in table The systematic erros have been determined 
using the same procedure as for the general fit. The correlation coefficients for this fit are listed in 
table The results are in agreement with the prediction of a V—A structure of the charged weak 
current. 


r —> ^ UiUr 

P 

0.781 ±0.028 ±0.018 

e 

0.98 ±0.22 ±0.10 


0.65 ±0.14 ±0.07 

0 

0.027 ± 0.055 ± 0.005 


Table 7: Result of the e-p, universality fit. 





V 

P 

-0.521 

-0.492 

0.348 



0.273 

- 0.022 

^5 



-0.094 


Table 8 : Correlation coefficients for the e-p universality fit. 
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4.3 Limits on the couplings 


From the measurement of the Michel parameters, limits on the absolute values of the couplings 
(see eq. |l]) can be extracted. This is done by constructing positive-semidefinite expressions from the 
measured parameters. A general approach to find such expressions is to use the boundaries of the 
physically allowed parameter space as shown in figure In three dimensions, the physically allowed 
region of the three parameters p, ^ and ^5 forms a tetrahedron p^ .P| 

An upper bound on |5 ^l| as well as weak upper bounds on I^^lI and I^^lI can be set from the 
expression 1 — p. Limits on and can be retrieved from the expression p — ^6. Limits on 

the remaining couplings |plrI I^l^rI follow from the probability for the decay of a right-handed r 
lepton which is given below. An even stronger limit on |plrI t)y regarding a plane in the 

3-dimensional parameter space {p,^,^S) which yields the expression 1 — p + ^ ^ ^ ^S. The explicit 

dependence of these expressions on the couplings is: 


^ — P = ^bLLl + I^LrI + I^RlI + ^Ifl'RRl + Y^I^LlI 
+ Y^bLR + ^PlrI + Y0 IS'RL + ^PrlI +Yglfl'RRl 
p-i^ = oI^RrI +oI5LR“2plr| ToblRl ) 




Is'lrI +yI^rrI + ' 


(26) 

(27) 

(28) 


Only the coupling |pllI cannot be constrained since it cannot be distinguished from the Standard 
Model coupling IpllI oR basis of the four Michel parameters. An upper bound on |pllI would require 


the measurement of correlations between one of the neutrinos and the charged lepton |47] or of the 
cross section for the process IHl • 

The probability that the r lepton decays as a right-handed particle can be calculated as the sum 
of all couplings of the type normalized according to equation 


l,.s ,2 , 1 , 


Qr — + Is'lrI + I^rrI +3|5'lr,| 


.V |2 


T |2 




(29) 


Using the correlations between the parameters the e-p universality fit yields: 


Qr = 0.089 ±0.131 < 0.304 (90%C.L.) 


(30) 


From the result of the global and of the e-p universality fit, respectively, the bounds given in table ^ 
can be set at the 90 % confidence level. Figure ^ shows the limits on the universal coupling constants 
normalized to their maximum values (^^ax = 2, ffmax = 1 ^^ax = ■^) introduced in section IT. 


4.4 Mass of a charged Higgs boson 

In models with two scalar field doublets, such as the Minimal Supersymmetric Standard Model 
(MSSM), the existence of a charged Higgs boson is assumed which contributes to the r decay through 
a scalar coupling. The value of the additional coupling is, assuming vanishing neutrino masses |48[-^: 


g f tan (5 

gr = -mi ITlr - 

VmH± 


( 31 ) 


^The three plots in figure H show the projections of this tetrahedron. 
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r — 

> e 

r —) 

pV^Vr 

r — 


^RR 

< 

1.36 

< 

1.25 

< 

1.05 

^LR 

< 

1.40 

< 

1.27 

< 

1.10 

fl'RL 

< 

2.00 

< 

2.00 

< 

2.00 

^LL 


< 2 


< 2 


< 2 

ffRR 

< 

0.68 

< 

0.62 

< 

0.53 

^LR 

< 

0.43 

< 

0.39 

< 

0.35 

fl'RL 

< 

0.56 

< 

0.55 

< 

0.52 

^LL 


< 1 


< 1 


< 1 

5lr 

< 

0.41 

< 

0.37 

< 

0.32 

S'rl 

< 

0.52 

< 

0.52 

< 

0.51 


Table 9: 90 % confidence limits on the coupling constants for the global fit and under the assumption 
of e-p, universality. As mentioned in the text, no limits can be set on the couplings gf^i^ and (7 ll which 
are listed only for completeness. 


Here m-H± is tbe mass of the charged Higgs boson, tan f3 is the ratio of the vacuum expectation values 
of the two Higgs doublets and ^ denotes either eoi p. Under the assumption that the neutrinos are still 
left-handed, the couplings are of the type After applying the normalization Ni = I^LLp+il^RR 
the Michel parameters can be written as:|^ 

3 

= 7 > 


= 


^ (5rr,£/2) 


1 + (5rr,£/2) 


(32) 




m = 


9RR,e/‘^ 


1 + (fi'RR,£/2) 


Using the above relations, the value of mjji / tan fi can be htted directly to the data. The likelihood 
function saturates for high Higgs boson masses or small values of tan (5. From the log likelihood a 
limit can be extracted as:[^ 


mf{± > 0.97 X tan/3 GeV (95% C.L.). 


(33) 


4.5 Left-right symmetric model 

In left-right symmetric models, parity violation of the charged current is caused by spontaneous 


symmetry breaking. In such models a second W boson is assumed The mass eigenstates 

Wi ^2 are not necessarily identical to the weak eigenstates Wl,r as mixing can occur. The model is 
parameterized by the mass ratio f3 of the physical eigenstates, 


/3 = ^ 
W 2 


m 


(34) 


'^^This can be verified by inserting the normalized ^ll and into the definition of the Michel parameters given in 
equation ^ 

^'^The following limits are determined as the value at which the log likelihood has dropped by the amount that 
corresponds to the quoted confidence level . 
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and by the mixing angle C that connects the physical masses to the masses of the weak eigenstates: 


"^Wl,2 




In this model the Michel parameters can be written as: 


p = - cos^ C ( 1 + C + 


4/3 


^ = COS^ ~ 0 


1 + /32 
1-/32 
1 +^’ 


tan^ 


{4 = jC 

p = 0. 


(35) 


(36) 


A limit on /3 can be transformed into a limit on mwj by using the direct measurement of the W mass: 
mwi = (80.43 ± 0.08) GeV @. 

For arbitrary mixing the upper plot of figure shows the one, two and three a contours of the log 
likelihood as a function of (5 and C,. Integration of the two-dimensional likelihood over C, yields the 
corresponding function shown in in the lower left plot of the same figure. From this likelihood a limit 
on mw 2 which is valid for arbitrary mixing can be extracted as: 

mw2>137GeV (95%C.L.). (37) 


Similarly, integration over mw 2 allows one to set bounds on the mixing angle independently of the 
W 2 mass (lower right plot): 

Id <0.12 (95%C.L.). (38) 

For C = 0, W 2 and Wr, become identical, and a limit on mwR can be given from a fit of the Michel 
parameter ^ alone. In this case there is no mixing but an additional coupling to a pure right-handed 
W that is proportional to its inverse mass:|^ 


V 

5ll,rr 




L,R 


(39) 


Under the assumption of no mixing the following limit is extracted from the log likelihood function: 


mwR>145GeV (95%C.L.). 


(40) 


5 Summary 

The Michel parameters of the leptonic r decays have been measured from the data collected with 
the OPAL detector in the years 1990 to 1995. The parameters pi, ^£, {^6)£ and (with i = e, p) 
were extracted from the energy spectra of the charged decay leptons and from their energy-energy 
correlations. A new method has been presented which involves a global likelihood fit of Monte Carlo 
generated events with radiative corrections at the generator level and complete detector simulation 
and background treatment. In the framework of the most general Lorentz structure for both leptonic 
decays the result of the global fit is: 

Inserting the normalized and ^rr into the definition of the Michel parameters (eq. yields the above relations 
for C = 0. 
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Pe = 0.779 ± 0.047 ± 0.029, 
= 1.13 ±0.39 ±0.14, 
(C5)e = 0.72 ±0.31 ±0.14, 


Pf, = 0.777 ±0.044 ±0.016 
= 0.79 ±0.41 ±0.05, 

= 0.63 ±0.23 ±0.09, 

7?^ = 0.010 ± 0.065 ± 0.001 


where the value of p is dominated by a constraint using the previously published values of the leptonic 
branching ratios and the r lifetime. The r polarization has been inferred from neutral current data. 
The Michel parameters have also been measured under the assumption of e-p universality and in terms 
of specific models. The e-p universilty fit yields; 


p = 0.781 ±0.028 ±0.018, ^ = 0.98 ± 0.22 ± 0.10, 

p = 0.027 ±0.055 ±0.005, = 0.65 ± 0.14 ± 0.07. 

Limits have been obtained on the individual coupling constants as well as on the masses of new gauge 
bosons, such as a right-handed W boson, Wr, and a charged Higgs boson. No indication for new 
physics processes has been observed. The results are in agreement with the V—A prediction of the 
Standard Model. 
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Figure 1: Selection efficiencies for the e-h (fi-h.) event samples as a function of the scaled lepton 
energy (momentum), xe (xp), as determined from the Monte Carlo after the preselection including all 
fiducial cuts. 
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Figure 2: The scaled r —> et'e^'r electron energy decay spectrum from e-h events. The quoted Michel 
parameters are the subset from the global fit that determines the plotted electron Monte Carlo spectrum. 
The dotted line represents the Standard Model expectation. 
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Figure 3; The scaled r —> 'muon momentum decay spectrum from /r-h events. The quoted 

Michel parameters are the subset from the global fit that determines the plotted muon Monte Carlo 
spectrum. The dotted line represents the Standard Model expectation. 
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Figure 4: The correlated r ^ ei^e’^r and r —> decay spectra from e-/i events together with 

the Monte Carlo distribution from the global fit. The spectrum of the scaled energy, xe, of the decay 
electron is shown in slices for each bin of the scaled muon momentum, Xp. The dotted line represents 
the Standard Model expectation. 
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Figure 5: The correlated r ^ and r —> decay spectra from e-/i events together with the 

Monte Carlo distribution from the global fit. The spectrum of the scaled momentum, Xp, of the decay 
muon is shown in slices of each bin of the scaled electron energy, xe- The dotted line represents the 
Standard Model expectation. (The plots show the same bins as figure Q m a different order.) 
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Figure 6: Presentation of the global fit result in 2-dimensional projections of the Michel parameter 
space. Shown are the 1 and 2 a ellipses of the fitted parameters. The outer (solid) line encloses the 
physically allowed region of the parameters. The dashed line confines the area that remains when the 
parameter which is not plotted is fixed to its Standard Model value. The closed circles indicate the 
basis spectra used in the fit. The labeled one represents the Standard Model. 
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Figure 7: 90% confidence limits on the normalized coupling constants g'luilgZiax (with = 2, 
ffmax = 1 ffmax ~ Under assumption of e-/i universality. The Standard Model coupling g^i^ 
which is not constrained is chosen to be real and positive. 
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Figure 8; Log likelihood as a function of the mass of a second W boson, and the mixing angle ( 

in the left-right symmetric model (upper plot). From this likelihood function limits on mw 2 (lower left 
plot) and ( (lower right plot) are calculated after integration of the likelihood over the second variable. 
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Type of interaction 
£-z/£-vertex (S) r-z^r-vertex 


Coupling constants 


Michel parameters 


p i V 
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0 

S-P g S+P 

9rr = 2 

3/4 

-1 
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p g p 

“sEl “ 9rl = 9lr “ “9rr = 1 

3/4 

0 

0 

0 

S+P g S+P 
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3/4 

-1 
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4l = 2 

3/4 
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T g T 
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II 

1/4 
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50% V-A g V-A, 50% S g S 

tT< 

tr 

CO 

"^1 
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50% V-A g V-A, 50% S-P g S+P 
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9ll = 9rr = 1/2, Srr = Sll = ^ 
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0 

0 
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50% V-A g V-A, 50% V+A g V-A 
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-1 
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50% V+A g V+A, 50% V-A g V+A 
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3/8 

1 

-3/8 
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50% V-A g V+A, 50% V+A g V-A 

9lr = 9rl = \/l/2 

0 

0 

0 

0 

67% V-A g V+A, 33% V+A g V-A 

9lr “ \/2/3, 5rl = \/l/3 

0 

1 

0 

0 

75% S+P g S+P, 25% T g T 

sEr = 4l = \/^- 9^ = 

1 

0 

0 

0 

12.5% S+P g S+P, 50% V+A g V+A, 37.5% T g T 

9lr “ 9rl “ 9lr = 9rl = 1/2, 9 ^ = 1/4 

0 

0 

0 
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Table 10: Example coupling constants and corresponding Michel parameters 


30 

















References 

[1] ALEPH, D. Buskulic et ai, Phys. Lett. B346, 379 (1995). 

[2] L3, M. Acciarri et al, Phys. Lett. B377, 313 (1996). 

[3] ARGUS, H. Albrecht et al, Phys. Lett. B246, 278 (1990). 

[4] ARGUS, H. Albrecht et al, Phys. Lett. B316, 608 (1993). 

[5] ARGUS, H. Albrecht et al, Phys. Lett. B341, 441 (1995). 

[6] ARGUS, H. Albrecht et al, Phys. Lett. B349, 576 (1995). 

[7] ARGUS, H. Albrecht et al, (1997), hep-ex/9711022. 

[8] SLD, K. Abe et al, Phys. Rev. Lett. 78, 4691 (1997). 

[9] GLEO, R. Ammar et al, Phys. Rev. Lett. 78, 4686 (1997). 

[10] GLEO, J. P. Alexander et al, Phys. Rev. D56, 5320 (1997). 

[11] Grystal Ball, H. Janssen et al, Phys. Lett. B228, 273 (1989). 

[12] W. Petscher, H. J. Gerber, and K. F. Johnson, Phys. Lett. B173, 102 (1986). 

[13] A. Pich and J. P. Silva, Phys. Rev. D52, 4006 (1995). 

[14] F. Scheck, Phys. Rept. 44, 187 (1978). 

[15] L. Michel, Proc. Phys. Soc. A63, 514 (1950). 

[16] G. Bouchiat and L. Michel, Phys. Rev. 106, 170 (1957). 

[17] G. A. Nelson, Phys. Rev. D40, 123 (1989). 

[18] W. Fetscher, Phys. Rev. D42, 1544 (1990). 

[19] D. Bardin et al, CERN Report No. TH 6443/92, 1992 (unpublished). 

[20] D. Abbaneo et al, GERN Report No. PPE-97-154, 1997 (unpublished). 

[21] OPAL, G. Alexander et al, Z. Phys. C72, 365 (1996). 

[22] OPAL, K. Ahmet et al, Nucl. Instrum. Meth. A305, 275 (1991). 

[23] OPAL, A. Honma, Nucl. Instrum. Meth. A348, 409 (1994). 

[24] OPAL, R. Akers et al, Z. Phys. C65, 1 (1995). 

[25] OPAL, G. Alexander et al, Phys. Lett. B266, 201 (1991). 

[26] OPAL, R. Akers et al, Z. Phys. C66, 31 (1995). 

[27] S. Jadach, B. F. L. Ward, and Z. W§s, Comput. Phys. Gommun. 66, 276 (1991). 

[28] S. Jadach, J. H. Kiihn, and Z. W§,s, Gomput. Phys. Gommun. 64, 275 (1991). 

[29] M. Schmidtler, Karlsruhe U., lEKP Report No. KA-93-14, 1993 (unpublished). 

[30] S. Jadach, M. Jezabek, J. H. Kiihn, and Z. Wcis, Comput. Phys. Gommun. 70, 69 (1992). 


31 



[31] R. Decker, S. Jadach, J. H. Kiihn, and Z. W§s, Comput. Phys. Commun. 76, 361 (1993). 

[32] E. Barberio, B. van Eijk, and Z. W§,s, Comput. Phys. Commun. 66, 115 (1991). 

[33] M. Bohm, A. Denner, and W. Hollik, Nucl. Phys. B304, 687 (1988). 

[34] F. A. Beerends, R. Kleiss, and W. Hollik, Nucl. Phys. B304, 712 (1988). 

[35] S. Jadach, B. F. L. Ward, and Z. W§s, Comput. Phys. Commun. 79, 503 (1994). 

[36] R. Battacharya, J. Smith, and G. Crammer, Phys. Rev. D15, 3267 (1977). 

[37] J. Smith, J. A. M. Vermaseren, and G. Grammer, Phys. Rev. D15, 3280 (1977). 

[38] OPAL, J. Allison et al., Nucl. Instrum. Meth. A317, 47 (1992). 

[39] R. Brun, F. Bruyant, M. Maire, A. C. McPherson, and P. Zanarini, CERN Report No. DD/EE/84- 
1, 1987 (unpublished). 

[40] F. James and M. Roos, Comput. Phys. Commun. 10, 343 (1975). 

[41] R. Barlow and C. Beeston, Comput. Phys. Commun. 77, 219 (1993). 

[42] A. Stahl, Phys. Lett. B324, 121 (1994). 

[43] W. J. Marciano and A. Sirlin, Phys. Rev. Lett. 61, 1815 (1988). 

[44] OPAL, G. Alexander et al, Phys. Lett. B374, 341 (1996). 

[45] OPAL, R. Akers et al, Z. Phys. C66, 543 (1995). 

[46] A. Rouge, Ecole Polytechnique Report No. X-LPNHE-95-01, 1995 (unpublished). 

[47] C. Jarlskog, Nucl. Phys. 75, 659 (1966). 

[48] H. E. Haber, SLAG Report No. SLAC-0343, 1989 (unpublished). 

[49] W. Hollik and T. Sack, Phys. Lett. B284, 427 (1992). 

[50] B. McWilliams and L.-F. Li, Nucl. Phys. B179, 62 (1981). 

[51] J. Polak and M. Zralek, Nucl. Phys. B363, 385 (1991). 

[52] J. Polak and M. Zralek, Phys. Rev. D46, 3871 (1992). 


32 



